M=0.1;            %nw/w

S=0:0.05:1;
kw=S.*S;
knw=2*(1-S)*M+(1-S).*(1-S)*(1-2*M);



% rel=[
% 3.544304e-01	5.855768e-01	4.150132e-01
% 4.873418e-01	7.393137e-01	2.608533e-01
% 6.424051e-01	8.734058e-01	1.261518e-01
% 7.143987e-01	9.180863e-01	8.117699e-02
% 1.645570e-01	3.037984e-01	6.971175e-01
% 8.860759e-02	1.704963e-01	8.304338e-01
% 3.876582e-02	7.643241e-02	9.244078e-01
% 2.294304e-02	4.559984e-02	9.551563e-01
% 2.500000e-01	4.396578e-01	5.611605e-01
% 3.037975e-01	5.176120e-01	4.830805e-01
% 8.860759e-01	9.770935e-01	2.240170e-02
% ];


% rel=[
% 3.544304e-01	2.028036e-01	4.164945e-01
% 4.873418e-01	3.138548e-01	2.625394e-01
% 6.424051e-01	4.683731e-01	1.280755e-01
% 7.215190e-01	5.598910e-01	7.859092e-02
% 1.645570e-01	7.767469e-02	6.982459e-01
% 8.939873e-02	3.892875e-02	8.299774e-01
% 3.876582e-02	1.613676e-02	9.255052e-01
% 2.294304e-02	9.408485e-03	9.562428e-01
% 2.500000e-01	1.291403e-01	5.624445e-01
% 3.037975e-01	1.654361e-01	4.844800e-01
% 8.109177e-01	5.952019e-01	1.197585e-01
% 8.908228e-01	7.417540e-01	4.795541e-02
% ];


%M=0.2
% rel=[
% 3.544304e-01	2.028050e-01	4.164893e-01
% 4.873418e-01	3.138551e-01	2.625394e-01
% 6.424051e-01	4.683822e-01	1.280757e-01
% 7.207278e-01	5.600643e-01	7.901070e-02
% 1.613924e-01	7.672413e-02	7.038990e-01
% 8.544304e-02	3.742969e-02	8.376744e-01
% 3.876582e-02	1.613678e-02	9.255059e-01
% 2.294304e-02	9.408474e-03	9.562429e-01
% 2.500000e-01	1.291412e-01	5.624236e-01
% 3.037975e-01	1.654359e-01	4.844777e-01
% 7.207278e-01	5.600643e-01	7.901070e-02
% 8.006329e-01	6.644903e-01	4.178846e-02
% ];



%M=0.1
% rel2=[
% 3.544304e-01	1.479560e-01	4.162578e-01
% 4.810127e-01	2.473024e-01	2.693216e-01
% 6.455696e-01	3.964427e-01	1.258444e-01
% 7.064873e-01	4.816169e-01	9.240379e-02
% 1.621835e-01	4.701621e-02	7.021438e-01
% 8.702532e-02	2.100401e-02	8.342377e-01
% 3.876582e-02	8.329947e-03	9.251491e-01
% 2.294304e-02	4.754495e-03	9.558007e-01
% 2.500000e-01	8.638723e-02	5.621657e-01
% 3.037975e-01	1.160116e-01	4.842648e-01
% 7.064873e-01	4.816169e-01	9.240379e-02
% 7.911392e-01	5.895518e-01	5.020975e-02
% %8.030063e-01	5.830017e-01	4.586742e-02
% ];
% 


% M=0.1;raidus 40
rel=[
1.612261e-02	3.306291e-03	9.682235e-01
2.906051e-02	6.201679e-03	9.429032e-01
6.329618e-02	1.488420e-02	8.773898e-01
1.046975e-01	2.767638e-02	8.016311e-01
1.580414e-01	4.761822e-02	7.090089e-01
2.243232e-01	7.794642e-02	6.013636e-01	
3.009554e-01	1.206412e-01	4.885154e-01
3.891322e-01	1.792956e-01	3.730189e-01
4.866640e-01	2.546374e-01	2.636316e-01
5.979299e-01	3.511920e-01	1.625381e-01
%7.221338e-01	4.765816e-01	7.961290e-02
%8.556927e-01	6.669537e-01	3.336153e-02
8.957006e-01	7.274614e-01	1.733430e-02
%8.564889e-01    6.714116e-01    2.341085e-02
8.207022e-01    6.786097e-01    2.723170e-02
7.136076e-01    5.204196e-01    1.189515e-01
8.989055e-01    6.874107e-01    1.071456e-02
];


rel=[0.0230708 0.00408246 0.958107
0.0644391 0.0143522 0.881101
0.117741 0.0367384 0.783984
0.198091 0.0507557 0.653255
0.297534 0.0982243 0.508841
0.416866 0.208452 0.339321
0.560859 0.33882 0.192312
0.631663 0.40974 0.135914
0.713604 0.502847 0.0816488
0.796341 0.604389 0.0415922
];




plot (S,kw,'r');
hold on
plot (S,knw,'g');
hold on
plot (1-rel(:,1),rel(:,3),'g^');
hold on
plot (1-rel(:,1),rel(:,2),'rx');
%================================
% hold on
% plot (1-rel2(:,1),rel2(:,3),'g*');
% hold on
% plot (1-rel2(:,1),rel2(:,2),'ro');
%================================
hold off
legend('WP analytical','NWP analytical','NWP LBM','WP LBM');


